Relaxation-to-creep transition of domain-wall motion in a two- 
dimensional random-field Ising model with an ac driving field 
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Abstract -Monte Carlo simulations of a two-dimensional, random-field Ising model with an ac 
driving field are used to study the relaxation-to-creep transition of domain-wall motion at low 
temperatures. The resultant complex susceptibility x = x' ~ exhibits features in agreement 
with the experiments of ultrathin ferromagnetic and ferroelectric films: The semicircle and straight 
line in the x'"X" V^ot are Cole-Cole signatures of relaxation and creep states, respectively. The 
exponent /? describing the creep motion is measured, and an intermediate state between the 
relaxation and creep states is identified. 



Domain-wall dynamics has become an important topic 
of study in magnetic devices, nanomaterials, thin films 
and semiconductors [IHl]. In particular, quenched ran- 
domness in ferroic materials, e.g. ultrathin ferromagnetic 
and ferroelectric films, fundamentally affects the response 
to an external field [5HTT| . From a purely theoretical point 
of view, the topic is also essential for understanding non- 
equilibrium dynamics in disordered media. For a dc (di- 
rect current) driving field, H , the domain- wall motion ex- 
hibits a depinning transition at zero temperature. The de- 
pinning field. Hp, separates the regimes of static pinning 
[H < Hp) and friction-limited viscous sHde {H > Hp). 
At low, but non-zero, temperatures, the sharp depinning 
transition is softened and a thermally activated creep state 
appears [HEllIlHig • 

Although the stationary states for a dc driving field arc 
now understood, theoretical results for the dynamic states 
in an ac(oscillating) driving field, H = Hq exp(i27r/t) is 
limited. At low temperatures, four dynamic states (relax- 
ation, creep, sliding and switching) may be classified in 
the so-called Cole-Cole diagrams [SI [71 [171 [18] . Segmental 
domain-wall relaxation without net wall motion occurs at 
high frequencies and weak driving fields [5l[7l [l"9H2T| . For 
low frequencies and strong fields, the domain- wall dynam- 
ics turns into periodic switching between differently poled 
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states [22]. Between them, two different dynamic states 
are complemented |20[|23j . Sliding of a domain wall is 
driven by an external field with Hq > Hp, while creep is 
the thermally activated motion for weaker fields, Hq < Hp. 
Experimentally, all four dynamic states are observed in 
ferromagnetic and ferroelectric materials [Tl fTTll^ . 

Theory suggested that a dynamic transition might sepa- 
rate the relaxation and creep states |13j . and very recent, 
experimental evidence was found in ultrathin ferromag- 
netic and ferroelectric films, as well as in liquid crystals 
[SllTllTS]. Theoretically, a Debye law and a power law may 
describe the complex susceptibility in the relaxation and 
creep states respectively [Tl lTOll^l^ . 

Current theoretical approaches to domain-wall dynam- 
ics in ultrathin ferromagnetic and ferroelectric films are 
typically based on the Edwards- Wilkinson equation with 
quenched disorder (QEW) [5l[2Ql[23l[25] . The QEW equa- 
tion is phenomenological, containing little microscopic in- 
formation, and describes neither the relaxation state nor 
the relaxation-to-creep transition [57] . To understand the 
domain- wall motion at a microscopic level, we need lattice 
models based on microscopic structures and interactions. 
The random-field Ising model with a driving field (DR- 
FIM) is a candidate [I2 1I28H3T] , at least to capture certain 
robust features of domain-wall motion, although it is only 
a minimalistic model and does not include all interactions 
in real materials. Very recently, the depinning transition 
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of the DRFIM model was examined with the short-time 
dynamic approach, and the critical exponents were ac- 
curately determined [16l|32]- The domain- wall roughen- 
ing process at the order-disorder phase transition has also 
been seen 

In this Letter, our goal is to simulate the relaxation state 
and relaxation-to-creep transition of the domain- wall mo- 
tion at low temperatures using the two-dimensional DR- 
FIM model. Numerical results will be compared with 
those experiments of disordered ferromagnetic and ferro- 
electric films. 

The two-dimensional DRFIM model is defined by the 
Hamiltonian 



(1) 



where Si = ±1 is a spin at site i of a square lattice, hi 
is the quenched random field uniformly distributed within 
an interval [—A, A], and H(t) is a homogeneous driving 
field. For our study, we typically take the coupling con- 
stant J = 1 and the disorder parameter A = 1.5 [16], 
and perform simulations at low temperatures, T, with a 
periodic field H{t) = Hq cxp{i2TTft). When Hq is suffi- 
ciently weak, the field strength dependence is negligible, 
so we usually fix Hq = 0.01. However, extra simulations 
are also performed with other values of the parameters J, 
A and Hq to examine their dynamic effects. The range of 
frequencies is lO"'' < / < 10 Hz. A total of 20000 runs 
were made so that statistical error bars in the results are 
usually smaller than the symbols used in the figures to 
follow. A total cpu time of 6 processor-years was used for 
the simulations. 

The initial state is semi-ordered with a domain wall in 
the y direction. Antiperiodic and periodic boundary con- 
ditions are used in x and y directions respectively. To 
eliminate the pinning effect irrelevant for disorder, we ro- 
tate the square lattice such that the initial domain wall is 
in the (11) direction [T2l[l6]. With this initial state, wc 
update randomly selected spins with the heat-bath algo- 
rithm. Denoting Lj. and Ly as the lattice sizes in x and y 
directions, a Monte Carlo time step is defined by x Ly 
single-spin flips. As time evolves, the domain wall moves 
and roughens, while the bulk remains almost unchanged 
since the temperature is low. The domain interface width 
grows slowly, so lattice size L^; = 12 and Lj, = 512 is used. 
Data for — 24 confirm that the finite-width effect is 
negligible. 

To describe the domain interface, we measure the mag- 
netization from the height function 
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Here Sxy denotes a spin at site {x,y), and < . . . > repre- 
sents both the statistic average and average in the y direc- 
tion. After the stationary magnetization hysteresis loop is 



obtained a.tt>to, the complex susceptibility x = x'~*x" 
is calculated by p51[?7] 



x{f,T) = 
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where P = 1/ f is the period of the ac driving field. In 
our simulations, to is about 20 periods, and another 20 
periods are performed for calculating x- To extract the 
creep motion induced by the temperature, we introduce a 
pure susceptibility 

Dx{f,T) = x{f,T)-x{f,T = 0). (4) 

For the weak driving field, M{t) oscillates by 



M{t) ^ Acos{2TTft - 5) + M. 



(5) 



Here A is the amplitude, 5 is the phase shift and Af is the 
average magnetization. Substituting Eq. ([5]) into Eq. ([3]), 
we find 



A{f,T)^2HoV. 



X'- 
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tan[<5(/,r)] = ^. (6) 
X 



Theoretically, the Debye law derived from non- 
interacting two-level systems may characterize the relax- 
ation state at T = or at high frequencies and T > 



Xoo(l//)^ 



X = , _2 ' \n 



(1//)^ 



where r is a relaxation time and Xoa = = oo)- 

It gives rise to a semicircle in the x'-x" (Cole-Cole) plot. 
With Eqs. © and ([7]), one may obtain 



A oc 



(1//) 



tan (5 oc (1//)" 



(8) 



At low frequencies and T > 0, an power law may describe 
the complex susceptibility in the creep state. If one as- 
sumes X = Xoo[l + (*/''') ~'^] with < /3 < 1, 



x'-XooOc(1//t)'5cos(/3V2) 
X"cx(l//T)'5sin(/3V2) 



(9) 



which is depicted by a straight line x" = (x' ~ 
Xoo)tan(/37r/2) in the Cole-Cole plot [5l[7]. In the low- 
frequency regime, xifjP = 0) w Xoo- Thus a power-law 
behavior of Dx is expected from Eq. The intermedi- 
ate regime between the relaxation and creep states looks 
complicated [Sl[71[251[171|3n] . 

In Fig. [U the time oscillation of the magnetization 
is displayed for Hq = 0.01, / = 0.1 and at T = 0.2 
in the upper panel. Compared with the external field 
H{t) = ffo exp(i27r/t), there exists a phase shift S as in- 
dicated in Eq. ([5]). Both the amplitude and phase shift 
of the magnetization are generally frequency-dependent, 
and could be also temperature-dependent. In the lower 
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panel, it is illustrated how the magnetization is driven by 
the external field. One can show that the area of the loop, 
S, is equal to x"- 

In Fig. [2] we show the Cole-Cole plot of the complex sus- 
ceptibility for the DRFIM model in a weak field Hq = 0.01 
at low temperatures. A perfect semicircle for the relax- 
ation state is observed at T = 0, consistent with Eq. (O 
and recent experiments [5l[7lfT8]. As T increases, the creep 
motion gradually dominates and a linear function is ob- 
served at T = 0.2 with a slope 0.76(3). The exponent 
13 = 0.41(1) is then calculated from tan(/37r/2) = 0.76(3). 
Two characteristic frequencies /i and /2 can be identified 
by the local maximum and minimum of x". In experi- 
ments, either /i or /2 is considered as the transition point 
[5][7]. In fact, these two frequencies separate three different 
regimes, the relaxation state / > A, creep state / < A, 
and intermediate state /2 < / < /i- 

In Fig. [31 x' and x" a-rc plotted as functions of 1//. 
Note that /2(T 0) = while fi(T 0) = 0.17 Hz. At 
T — 0, the Debye law for the relaxation state in Eq. ^ 
is shown with solid lines. The relaxation time r « 7 is 
obtained from the fitting, equal to l//i. More accurately, 
a logarithmic correction is detected for x", which has been 
reported as a non-Debye-type contribution [25j . In the re- 
laxation regime, the curves of different T almost overlap, 
indicating that the relaxation motion is T-independent. 
In the creep regime, power-law behaviors of x' and x" in 
Eq. ^ are evident; however, measurements of the expo- 
nent /? are not precise, partly due to the existence of the 
relaxation motion at T = 0. Additionally, the exponent 
(3 estimated from x" is smaller than that from x' and in 
Fig. El 

To identify the pure creep motion induced by the tem- 
perature, we plot Dx' and Dx" as functions of 1// at 
different temperatures in Fig. |4l Power-law behavior 
is observed, and the slopes of Dx' yield /? = 0.23(1), 
0.23(1), 0.26(2) and 0.38(1) for T = 0.025, 0.05, 0, 1 and 
0.2, respectively. /3 = 0.38(1) at T = 0.2 is consistent 
with 0.41(1) obtained in Fig. [21 A power-law correction 
may extend the fitting to the numerical data. Similarly, 
from Dx", we find /3 0.013(2), 0.027(4), 0.106(11) and 
0.288(12), respectively. The difference of the exponent /3 
measured from Dx' and Dx" modifies the empirical ansatz 
in Eq. ®. ^(T ^ 0) « 0.2 of Dx' and /3(T 0) « 
of Dx" indicate that the creep motion induced by the 
temperature is different in x' and x" ■ Finally, an ap- 
proximate relation Dx{T) oc T is detected in the creep 
regime when T is low enough. According to Eq. ([9]) , this is 
qualitatively consistent with the temperature dependence 
of the relaxation time reported in the literatures [IHIISS] > 
T ~ exp{c/T) . 

To capture the features of the intermediate state be- 
tween the relaxation and creep states, we examine the dy- 
namic behavior of the amplitude A and the phase shift 
tan 5 of the magnetization, as shown in Fig. [51 At T = 0, 
Eq. ([8]) derived from the Debye law fits the data, and in- 
cluding a logarithmic correction improves the fitting of the 



phase shift tan (5. At T = 0.2, the tail of tan (5 in the creep 
regime exhibits power-law behavior with a nonzero slope 

0. 12(1), which reflects the difference of the exponent /3 
for x' and x" obtained in Fig. [31 Interestingly, we observe 
that the frequency /i characterizes the transition behavior 
of A, while A depicts that of tan S, i.e. the relaxation-to- 
creep transition occurs in two stages, somewhat similar to 
the scenario of the two-dimensional melting [36] , where an 
intermediate state exists between the solid and liquid. In 
the relaxation regime, i.e., / > /i, the domain wall re- 
laxes microscopically without macroscopic movement, the 
amplitude A is T-independent, and the phase shift tan 5 
decreases with 1// although it is weakly T-dependent. For 
/ < /i , the creep motion induces the T-dependent ampli- 
tude A. Within the intermediate regime, i.e., /2 < / < /i, 
however, the phase shift tan i5 decreases continuously with 
1//, and the domain wall still looks partly like the relax- 
ation state, especially at low temperatures. At / = /2, the 
phase shift tan S reaches its minimum. In the creep regime, 

1. e., / < /2, the domain wall creeps macroscopically, the 
amplitude A shows a linear behavior A{T) « aT -I- 5 at 
low temperatures, and the phase shift tan (5 has a power 
law dependence on 1//. 

In the literatures, the average magnetization M = 
{01/271) f M{t)dt is often referred as the dynamic order 
parameter [7l[22|, and the area of the magnetization hys- 
teresis loop S = (a;/27r) / M{H)dH is also concerned [18] . 
By the definitions, S ^ x" can be deduced, while M and 
A show similar dynamic behaviors. From our numerical 
simulations, the amplitude A and phase shift tan 5 may 
characterize the intermediate state more clearly. 

In the creep state, it is well know that the relaxation 
time is induced by the metastable state of the domain 
wall, i.e., T ~ exp{EB /T) with Eb being the energy bar- 
rier [26] . For the relaxation state, the Debye law in Eq. ([7]) 
can be derived with a simple two- level system [26| . The 
only time scale t in this case is still assumed to be the 
relaxation time of the metastable state. Such an assump- 
tion should be reasonable at not very low temperatures, 
and is supported by experiments [5l[37]. At zero or very 
low temperature, the scenario can be different, at least 
from the view of the DRFIM model. At zero tempera- 
ture, for example, we detect a relatively small relaxation 
time as shown in Fig. [3] The domain wall relaxes and 
oscillates between the stationary states corresponding to 
the dc driving fields H = ±ho and H — 0, and the re- 
laxation time is not induced by the metastable state with 
the energy barrier Eb- Qualitatively, it may be estimated 
by suddenly quenching the domain wall from the state of 
H = to H = ho, and measuring the relaxation time 
of the domain-wall velocity. From our numerical simula- 
tions, we observe that the relaxation time r increases with 
ho, while decreases with the coupling constant J and the 
disorder parameter A. Therefore, the temperature depen- 
dence of T varies in the parameter space oi ho, J and A. 
To well describe the delicate behavior in the experiments 
[5l[37], however, more realistic interactions should be in- 
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eluded. 

In summary, Monte Carlo simulations elucidate domain- 
wall motion in the two-dimensional random-field Ising 
model in an ac driving field. The relaxation and creep 
states are bounded by difi'erent characteristic frequencies 
/i and /2 of the ac driving field: At / = /i, the am- 
plitude A transits from T-independent to T-dependent, 
while at / = /2 , the phase shift tan S reaches its minimum. 
The complex susceptibility x = x' ~ *x"j the amplitude 
A and phase shift tan (5 of the oscillating magnetization 
were examined at low temperatures. The Cole-Cole plot 
agrees characteristically with experiment, and power-law 
behaviors of Dx'i^/f) and Dx"{^/ f) in the creep regime 
were observed. At T = 0.2, the exponent (3 = 0.38(1) 
obtained from Dx' is consistent with /3 = 0.41(1) mea- 
sured from tan(/37r/2) = 0.76(3); however, /3 = 0.29(1) 
estimated from Dx" is smaller and modifies the empirical 
ansatz in Eq. (|9]). Our values of the exponent (3 can be 
compared with those from experiments, e.g., j3 = 0.35(2) 
and 0.4 in Refs. 0110]. 
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Fig. 1: The upper panel shows the time oscillation of the mag- 
netization (circles), in comparison to the external field (the 
sohd hne) with Ho = 0.01, / = 0.1 and at T = 0.2. The 
lower panel demonstrates how the magnetization is driven by 
the external field. The area of the loop, 5, is equal to x" ■ 




Fig. 3: The spectra of x s-nd x" on a log-log scale. Solid 
lines show the Debye law and its logarithmic correction for the 
relaxation state, and the dashed line indicates a power-law fit. 
Frequencies /i and /2 are marked for T = 0.2. 




Fig. 2: Cole-Cole plot of the complex susceptibility for different 
temperatures and frequencies of the ac driving field. The solid 
line shows a semicircular fit, and the dashed line shows a linear 
fit. The characteristic frequencies/i and /2, marked with up 
and down arrows, denote the local maximum and minimum of 
X". Pluses show data for = 24 at T = 0.2. 




Fig. 4: The spectra of Dx on a log-log scale. Dashed lines 
represent power-law fits, and solid curves include corrections 
y = ax^^{l — c/x). The inset shows the spectra of Dx" ■ 
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Fig. 5: The amplitude A and phase shift tan 5 for the magne- 
tization as functions of 1/f on a log-log scale. The dashed line 
is a power-law fit, and the solid lines show the Debye law and 
its logarithmic correction in the relaxation state, /i and /2 are 
the transition points for A and tan S respectively. 
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